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ABSTRACT 


f 


Existing finite element programs PLISOP and PLIMEG are 
modified and adapted to the successful determination of 
stress concentration factors for thin elastic bars containing 
V-shaped round bottom notches and subjected to tensile or 
bending loads. Accuracy and reliability are demonstrated. 
Comparisons are made with previous analytical and experimental 
Brrerminations. Ihe effects of finite bar thickness are 


discussed. 
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СЕРЕ PRODUC TION 


Stress intensifiers, in the form of holes, grooves, key- 
КОО threads, ete., result in local modification of simple 
stress patterns and give rise to stress concentrations. The 
Gesigner or analyst is interested in determining the locally 
maximum, governing or limiting value of stress, designated by 


the symbol OF пае τς Calculated as the product of a 


ax 


шошу паї stress, c I CIO Op C Onceoncpation factor, 


nom? 


designated by the abbreviation SCF, thus: 


σ 3 SUE x O 
n 


max om 


An elementary formula is used to compute Gone the terms 
of which must be precisely defined; this accounts for the 
Bross geometry of the member and the loading which is applied. 
The SCF is intended to provide a corrective multiplying 
factor; its value depends on the details of local geometry 
Zur usually, on the particular type of loading that is 
applied. 

The designer employs available tabulations of SCF's, 
together with the appropriate and easily calculated om? 
Ploy ing the formula given above, to estimate the maximum 
n SS Which he compares with an allowable value; the latter 
Meme established by experience, by code, or by contract. 


The task of the person who supplies values of SCF is to 


determine values of c and of G ang there quotient; the 
max non 





SCF's. To make the values of SCF most widely useful, the 
Stress analyst selects an appropriate set of parameters to 
describe the geometry of a class of stress intensification 
problems. A number of individual problems are undertaken, 
taking care to vary the parameters so as to cover the useful 
range of each and thus to permit making reasonable inter- 
polations for intermediate values. 

The SCF's obtained in this manner are usually presented 
in the form of parametric curves drawn through the actual 
data points established. The designer then must interpolate 
between the curves to match the actual parameters of his 
Problem. 

The particular classes of problems selected for this 
UE), because of their prominence and importance, are 
characterized by a thin, flat bar under tension and/or 
mending. We consider two forms of stress intensifiers: 

ieee lwo notches On Opposite edges of the bar, located 
sxuEercecM about the long axis of the bar and 


shaped symmetrically with respect to a plane normal 
Eo Clie long axis of the bar. 


ОО А notch on only one edge of the bar. 
ШОО ОСОО cheonotehes are circular at the bottom. Details 
of the notehes, with sketches, are presented in Section III. 
ΜΝ. οἱ what follows, however, is applicable to other 
meometries. 

We confine attention to the two-dimensional case, i.e., 
the case of plane stress. Realizing that a condition of 
pone stress is never actually found in a real world 
МЕ састоп, the cfrect of three dimensionality on SCF's 


was also investigated, see Appendix C. 
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IA DAS NOME TAODS OF OBTAINING 5СЕ'5 


A. CLASSICAL ANALYSIS 

Neuber [11]has made the most extensive analytical 
investigation of SCF's. His book contains treatments of a 
number of geometries and loadings and his methods include 
a few cases of exact analytical solutions together with a 
ШШШ Бег of ingenious approximations and interpolations. The 
best known of Neuber's analyses is that of the deep hyperbolic 
ouch. 

1. Neuber's Notch Analysis 

Ὁ πο silinoarrfeoordinates and complex stress 

functions, Neuber presented a solution for axial tensile 
loading (without bending) of a geometry consisting of 
symmetrically located, hyperbolic shaped notches in an 


Aute sheet of unit thickness (Fig. 1). 





Da Ww EsSwuaetrlcehyperbolic notches in an infinite 
Sheer. 


an 





He arrived at the following solution for the stress 


Boncentration factor for this notch: 


4 cot n 
SCF = ο 0 
π - No т sin(2ng) 


where No η ΠΕΤΕΡ the particular notch 
Beometry, namely, the angle between asymptote anar novena 
centerline. The nominal stress is given by the axial force 
divided by the minimum width of the sheet at the notch. The 
hyperbolic boundary is described by the equation: 


E. x _ 


2 2 пра 
a cos np a sin ng 


1 
ll 
ET 


Ἱ halt’ width ac the bottom of the notch is р, The 


ШЕ тит radius of curvature is 
Е 2 
DEL SPEO (ng) 


lhe width of the sheet approaches infinity on each side, 
mm ke cases Of practical interest in which the notch is an 
Eucentation in a bar of finite width. Accordingly, for 
Neuber's deep hyperbolic notch, it is not possible to specify 
a Enoss заса ог the ratio of gross width to net width as 
is usually done for practical cases. 

Neuber further presents an approximation for shallow 


Eccches 3n a bar of finite width with notch depth specified. 


ke 





The analysis involves employing a solution for an elliptic 
mole in a bar. Somewhat arbitrarily the bar is 
cut through the center of the hole parallel to the edges of 


the bar and reassembled to form the notched bar (Fig. 2). 





Fig. 2. Neuber's shallow notch approximation 


The boundaries of the notched bar are not unloaded as they 
Exod be and there are unaccounted for dislocations along 
ШО axial, center line. 

Taking the half width of the bar as W and using D, 
Band © om as defined for the hyperbolic notch analysis, 


Neuber arrived at the following approximation for stress 


concentration factor for the shallow notch: 
12 
ο το 2 (W-D)/p]° 


Noting that nominal stress was defined the same way 


ШО cases, 1.e., on the basis of net area, Neuber devised 
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a nomographic scheme to obtain SCF for a wide range of the 


parameters: 
p = minimum radius of curvature 
Pe=ehelt width ef the bar at the notch (net half width) 
ү = half width of the undisturbed bar (gross half width) 


The nominal stress is с Е =(axial force per unit thickness)/2D, 

A summary of Neuber's solution for the hyperbolic 
notch is shown in Appendix B. 

2. Limitations of Neuber's Analysis 

The hyperbolic notch analyzed has a definite shape, 
namely a hyperbola having the asymptotes y = * x tan No ° The 
values presented by Neuber are based on two extremes and an 
interpolation between them. It must be recognized that other 
geometries only approximate these notch shapes and therefore 
Mies data must be judiciously applied. 

The reasoning used to arrive at this solution is 
ingenious but experimental evidence described later indicates 
a systematic difference between Neuber's analysis and the 
L vor of actual test sections. 

ERE Uo) арту of Neuber's Solution 

A designer may use Neuber's equations and perform the 
ШО сатгеа interpolation to determine a particular SCF based 
Oy the Parameters of interest. However, much of that work 
has already been done. 

The most extensive, readily available tabulation of 
SCF for this case is presented by Peterson[12]. He presents 


a family of curves using dimensionless ratios of the 


14 





parameters as arguments. Other authors, such as Faires [6] 


present the same data over a narrower range of parameters. 


B. EXPERIMENTAL INVESTIGATION 
Values of maximum stress may be determined by experimental 
stress analysis including strain gage measurements and photo- 
Stic evaluations. 
1. Strain Gage Measurements 
Kikukawa [9] used resistance strain gages to study 
Mies stress in a bar with two symmetrically located U-shaped 


merenes Subjected to an axial tensile load. 





T 


Fig. 3. Symmetric U-shaped notches 





The geometry of the test section (Fig. 3) is specified by: 


I 
1 


F ST ler bottom of the notch 


J 
II 


net half width 
W = gross half width 
Only the gross width of the bar was varied during the tests. 
Peeve lois 0: Strain Gage Measurements 

In using strain gages, difficulty arises from the 
ENehliy localized character of “the stress distribution at the 
place in question. A very small gage length Is necessary 
K cI satisfactory results. Accordingly, the physical 
Size of available strain gages limits their usefulness to 
the investigation of relatively large notches. Furthermore, 
models must be very carefully machined to obtain dimensions 
within satisfactory tolerances. Good adhesion of the gage to 
the model and correct placement and orientation on the specimen 
meee cLificult to obtain. 

Data scatter dictates that Some sort of smoothing 
process must be used before the results are useful to a 
designer. Large expenditures of time and other valuable 
resources are required to obtain even a few benchmarks which 
can reasonably be used to infer SCF's over some useful range 
Of the parameters employed. 

3. Photoelastic Evaluations 

Flynn and Roll [7] used photoelastic techniques to 
ο παν the stress in a bar subjected ze an axial tensile load 
for both U-shaped and hyperbolie-shaped symmetrically located 
notches. The geometry of the U-notched test section (Fig. 3) 


moe specified by: 
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“ο 
it 


radius at the bottom of the notch 


oO 
l! 


net half width 
W 


gross half width 





Fig. 4. Symmetric hyperbolic notches 


The geometry of the hyperbolic-notched test section (Fig. 4) 


is specified by: 


р = minimum radius of curvature 
D = net half width 
W = gross half width 


The gross width of the bar was the only parameter which was 
varied in these tests, also. 

Appl and Koerner [2] used photoelastic techniques to 
study the stress in a notched bar subjected to an axial load. 
whe notches were symmetrically located, V-shaped with various 


meles and had round bottoms. 


Es 








hic ὁπποινιο V Shaped notches 


The geometry of the test section (Fig. 5) is specified by: 


ti 


mace ate the bottom of the notch 


mete пате width 


I 
lI 


gross half width 

8 = half included angle 
Mey present SCF's for several notch angles using the ratio 
D/W as a parameter. 

Cole and Brown [5] used photoelastic techniques to 
study the stress in a bar with a U-shaped notch on OnLy “Orme 
side of the bar. An axial load was applied such that its 
line of application was through the center of the net cross 


section. 
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o 
J 1 < 





Fig. 6. U-shaped notches, one side only 


The geometry of the test section (Fig. 6) is specified by: 


R = radius at the bottom of the notch 
= net width 
W= gross width 


Mey present SCF's for several combinations of the parameters. 
uy Limitations of Photoelastic Evauluation 
Models, no matter how thin, do have a third dimension. 
These experimenters have chosen thin models purposely so that 
their results could Dc reported as plane stress values. 
Appendix C discusses the effect of the third dimension, showing 
that the stress at the base of the notch actually varies through 


tbe thickness. Photoelastic evaluation reveals the stress 
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averaged through the model thickness. The difference between 
the averaged and true maximum is small, but the three 
dimensional effect should be recognized. 

Errors due to several edge effects including those 
caused by residual machining stresses and time dependent 
material changes must be recognized. 

The same limitations regarding model preparation, 
Gata scatter and resource expenditure tnat were mentioned as 
limitations to strain gage measurements apply also to 


photoelastic evaluations. 


C. NUMERICAL TREATMENT OF ELASTIC FIELD EQUATIONS 

Several numerical techniques, notably the Finite Difference 
NOG, are available and have been used to solve partial 
differential equations. However, useful results to the 
elastic field equations are difficult to obtain, especially 
if the geometry of the problem is not very Simple. 

Recently a new numerical method called the Finite Element 
Method (FEM) has been extensively developed and used. In many 
ways it is similar to earlier numerical methods in that the 
partial differential equations are repiaced by a finite system 
of algebraic equations. There is, however, a big difference 
in the detail of this replacement and this difference in 
detail permits the user to deal with complex geometric 
shapes. 

The method has been well proved by the solution of 


problems which have elassically known solutions. Thus it 





seems to provide a new and powerful way to determine stress 
concentratlon factors. 
Considering all of the above, the Finite Element Method 


(FEM) was chosen for use in this study. 


al 





i DT ROBLEM SPECIFICATION 


A. GENERAL 

Before proceeding with FEM program selection, it will be 
convenient to have the problem specified in full detail. 
This is done So that further reference to computer program 
selection and modification can be with the requirements of 
mas problem specifically established. 

ШОО СООП £ we limited the Study to a family of thin, 
Mao bars with general end loading. The stress intensifiers 
were limited to two forms: a single noteh aná two symmetrically 
located notches. Now specifically we deal with that general 
ruine and a detailed description of notch Shape. A single 
meee shape is considered: the bottom of the notch is 
Circular, characterized by its radius, and the sides are 
V-shaped. The included angle of the V may be between 0° 
(U-shaped notch) and 180° (unnotched bar). 

This problem family was chosen because of its practical 
importance and because it includes much of the available 


experimental data. 


КОО FOUR PROBLEMS 

Division of the problem family into two general classes 
Based On the number of notches has already been Μπο σα σα. 
The two classes are shown in Fig. 7. 

We speak of general loading at the ends of the bar 


МОЕ ут the stress intensifier. However, by virtue of the 
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Princeıiplesorsanntaienant, if the bar is several times longer 
ΠΠ ICS IOC, as we assume and as is the case in practical 
Brolıleations, then only the static properties of the load, 
1.6. the net axial force and the moment (about, say, the 
Mentral axis of the bar) are of importance. Thus what we 
have calleá a general loading can be treated by the super- 
"ὅσα. in proper proportion, of a case with an axial force 
and no moment and a caSe with moment and ο τοι axial load. 
Based on this, the four specific problems shown in Fig. 8 


and described in Table I were selected for study. 
TABLE I 


που) сс быа оа 


Problem Number Notches 
[oa s 
2 loment 
3 


i 







Loading 









Axial Force 






The general loading of the class of problems with two 
notches can be obtained by an appropriate superposition of 


problems 1 and 2. The maximum stress (σα E) always occurs 


a 
C i Dase of one of the two notches and since they are 
ИСЕТ it 15 apparent that the maximum stress can always 


be calculated. 
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The general loading of the Glass Of problems with one 
notch can similarly always be obtained by an appropriate 
superposition of problems 3 and 4. However the maximum 
Stress does not always occur at the base of the notch. 
Consider first the bending of problem 4. Maximum tensile 
Stress Occurs on the straight boundary opposite the notch and 
maximum compressive stress at the bottom of the notch. 
Superimposing problem 3 can bring, the stress at the bottom 
{6 notch to any value, say zero, producing maximum stress 
em che boundary opposite the notch. 

Accordingly, for this class, data should be generated and 
made available both for the point at the bottom of the notch 


and also for the opposite point on the straight boundary. 


C. ONE STRUCTURE 

By taking advantage of symmetry, the four problems may 
be reduced to four cases each of which is completely described 
Buche single geometry of Fig. 9 and one of four different 
Bor boundary conditions. Face B is a circular aro of 
Eus Е. Face C is a straight line inclined at angle 6. 
Mime Other four faces, A, D, E and F, are straight lines 
parallel to the x or y axis. The length of face A is 
designeted D, the length of face E is designated W, and the 
length of face F is designated L. We take L = 5W so as to 
meee, Dy 5t. Venant's principle, that fine details of 
ШОШ саа пе do not influence the stress distribution in the 


Mmemeened portion of the bar. 
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y 1 


Face C 





Fig. Structure faces 


The geometry of the structure can now be adequately described 


by the dimensionless ratios Ü 5 о апа ө 


р. FOUR SETS OF BOUNDARY CONDITIONS 
The four cases to be described here and summarized in 
Tables IT and III are equivalent to the four EMT 
deScribed earlier. Cases c and d are problems 3 and ! 
directly. Appropriate superposing of two cases a and of 
two cases b along face F yield problems 1 and 2 respectively. 
ος pOMmt We delineate the boundary conditions 
corresponding to each of che, ους cases. Ie iS evident thar 
we have to specify values of the displacements u and v along 


КОО A. B, C, D, E and E and/or the values of the forces 


za 





N and ὃν applied along these faces by external agencies. 
Along faces B, C and D there are no externally applied forces 
and the displacements u and v are not specified (that is, 

they are unknowns which will be solved for by the FEM 
program). All internal points of the structure, likewise, 
have no specified values of u or v and are not externally 
loaded. 

Because of the symmetry in all four cases, we specify 
that u = 0 all along face A. Along face F, for cases c and 
there are no boundary forces and values of u and v are 
ШО О е Tor case a, face F is along a line of symmetry 
and we take v = O all along face FP. For case b, the condition 
SI antisymmetry is satisfied by taking u = O at every point 
along face F. In neither of these ΤΕ any 
ppeseribed boundary forces along face F. In all four cases, 
the point which marks the intersection of faces A and F is 
restrained so that u = v = 0. 

Pedane oare applied on face E. The displacements along 
face E are unspecified. All forces FY are zero. Values of 
ο along this face are such that they produce the required 
loadings as shown in Table II. 

A summary of boundary conditions for faces A, B, C, D 
End P is given in Table III. The conditions are given as 
triples (ανν Ες). All values of F. are zero. Ап ο. (*) 


y 
indicates that no value is specified. 
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TABLE II 


EXTERNAL LOADING ALONG FACE E 


Static Effect 
Case 





The displacements along face E are unspecified. All forces 
P are zero. 
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TABLE III 


BOUNDARY CONDITION SPECIFICATION 
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IV. SOFTWARE SELECTION 


A. BASIC PROGRAM SELECTION 

The purpose of the study was not only to obtain some useful 
data (i.e., SCF's) to supplement and refine that already 
E uable through other sources, but also to perfect computer 
Ой аге Гог generating Such data and to establish its 
reliability. Reliable, efficient programs capable of solving 
the two-dimensional plane stress problem and three-dimensional 
Stress problem have been written and are available. Some of 
them are proprietary but many appear in the literature and 
οπου OL Others have been published. Unfortunately, mainly 
due to the state of the art in the field of computer science, 
Is (programs are machine dependent and therefore could not 
F considered for use in this investigation. 

The two programs selected were developed at the Naval 
Postgraduate School. PLISOP [4] is a two-dimensional finite 
element program capable of solving the plane stress problem. 
IT Utilizes isoparametric elements which means that curved 
boundaries can be represented. All data handling in this 
ram is contined to the main core storage of the computer. 
The program has internal documentation which includes instruc- 
tions for its use. TRISOP [0] is a three-dimensional finite 
element program utilizing isoparametric elements. Because 
of the large amount of data implicit in the three-dimensional 


E cUvsis. much of the program coding deals exclusively with 


E 





Шела шр the dala Most of the storage is outside the 
computer main core on direct-access devices which provides for 
the storage of large volumes of data. The price paid for 
this expanded storage is the large computational time needed. 
Шо шл documentation includes instructions for its use, 
Mor LO Llie start Of this investigation both of these 
programs had been debugged on the IBM-360 computer available 
"προ τη this sttay arfa™~solttions “to a number οἱ «problems 
had been compared to the classically obtained solutions to 
demonstrate that the programs are capable of producing very 
good approximations when used correctly. The remainder of 
this section deals specifically with the two-dimensional 
program. Cf. Appendix C for remarks concerning the use of 


W TSOP. 


B. MESH GENERATION 

The first step in employing a FEM program for the solution 
era Specific problem is to provide appropriate input data, 
eescribing the problem, in the form called for by the 
specifications of the particular FEM program employed. This 
requires the preparation of a very large volume of data. 
Material properties and program options must be selected, 
Of course, but the great preponderance of the input has to 
ο Πιτ geometry. The user must discretize the structure 
Dor elements and describe these elements and their inter- 
Eennection to the program. Coordinates for each node used 
(specific points on the element boundaries) must be provided 
along with load boundary conditions. The displacement boundary 


Меса Ons comstitute the third large block of information 


SE 





required. The elements and nodes are numbered for identi- 
mes tion using numbering conventions dictated by the coding 
of the FEM program. Specified data card format must be 
adhered to rigidly. The discretization ultimately selected 
for most of this study would require over 600 data input 
cards, 

Manual preparation of such a volume of data is unthinkable. 
However this problem is not new and others have dealt with it 
before and have provided programs to accomplish the task. 

Such a program is called a mesh generator. It is written to 
eceept a gross description of the structure under consideration 
(Fig. 12) and its basic function is to superimpose a grid, 
establishing the elements (Fig. 11) required. Coordinates 
fee! nodes are calculated. Elements and nodes are numbered 
according to the criteria of the FEM program which the mesh 
generator serves. The program output is in the form of a 

Beck of punched data cards. The deck, compatible with the 

FEM program, is in two parts: element descriptions (element 
connectivity cards) and nodal coordinates (element node cards). 

μα ασ а computer program originally written by J.R. 
Adamek as part of his NPS thesis [1] for the purpose of 
EEuIifving the preparation of input to PLISOP. Modifications 
to Adamek's program vere made as a part of the present study 
in order to make it more useful for the MM ode O1 πο ισα 
data for the notch stress intensification problem and making 


ΠΠ use of the computer facilities in solving the problem. 
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The basic input to PLIMEG is a description of the geometry 
оре problem, divided, by the author, into several large 
ENDOIvTISTons, referred to as superelements, or, briefly, as 
EE S. "The coordinates of the SEL corner points are specified, 
Em connectivity is indicated, and, in an appropriate fashion, 
the geometry between corner points is specified. Other inputs 
LO PLIMEG are m and n, the number of subdivisions required 
along two adjacent SEL boundaries. The numbers m and n must 
memconsistent on any boundary common to two SEL's. 

Internally, the program PLIMEG now establishes additional 
mega! points within each SEL, dividing each SEL intom xn 
smaller elements — the finite elements themselves. The output 
EU deck is compatible with PLISOP. A program option may 
ΙΝ. CG obtain 2 graphical output depicting the completed 
id. 

The subdividing done in PLIMEG is actually accomplished 
in a non-dimensional coordinate system (E,n) and the grid 
established there is mapped into the SEL by a coordinate 


transformation. 


See MODIFICATIONS TO PLIMEG 
i Stress Gradient Considered 
The structure considered in the study (Fig. 9) is of 
πι... εεαιδε οἱ the stress concentration. That concentra- 
n implies a large gradient which, in turn, demands small 
meements in the finite element discretization. The use of a 
MeeetOormly fine mesh is undesirable from the standpoint of the 


Ember of equations that must ultimately be solved. Not only 
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is the solution more time consuming butthe maximum number 

of elements is directly dependent on machine storage capacity. 
The IBM-360, using PLISOP, imposes an upper limit of slightly 
over 180 quadratic elements, and even that size of problem En 
only be considered by special treatment. 

Noting; that the high stress gradient is very localized, 
it is possible however, to use small elements in that 
Mimediate area and larger elements where the gradient is 
smaller. By extension, bands of different size elements 
may be used. This technique has been successfully employed 
by many users. Pfeifer [13], for example, uses the three- 
dimensional stress analysis program (TRISOP) in this way to 
Сее the problem of a concentrated force on the boundary of 
a semi-infinite body. His results compare well with Boussinesq's 
ENSssical solution. 

Accepting the usefulness of graded element size, a 
logical next step is to devise a way to vary element size 
progressively as the stress gradient changes. This vas 
Eompiished by modifying PLIMEG. Entering the program at 
the point where the nondimensional grid has been established, 
the element corner nodes, represented by (Ey n4) can be 
systematically moved by use of the relationships: 

expla; (éj -Enin?] 7 1 
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ma man 
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A wide range of element size variation is provided for by 
allowing the user to Specify the exponential coefficients 
а; апа an as PLIMEG input values. The intermediate nodes are 
M erted on the new grid lines thus established, equally 
spaced between corner nodes (one if quadratic elements are 
used and two if cubic elements are used). From this point 
mee cOOrdinate transformation and other details of PLIMEG 
are as originally programmed. Mesh generation as originally 
Procrammed and as modified is illustrated in Fig. 10 for 
a representative general SEL. Note that PLIMEG, as modified, 
meovides the capability of dividing a SEL into any specified 
number of increments with any specified exponential coefficient 
U I live, negative, or zero) in both the & and n directions, 
each completely independent from the other. 

Eee curved Boundary Node Relocation 

Isoparametric elements are used both for the SEL 

description and the elements in the completed mesh. Three 
element types are used. Each is capable of representing a 
different degree of variation along an element boundary. 
The element types, with their common names and nodal description 
are shown in Table IV. The maximum degree of polynomial which 
Beer exactly represented by each element is given. If 
straight line, parabolic, and/or cubic boundaries determine 
the geometric shape óf the structure, they may be exactly 
Eu c-nteg by proper choice of SEL type. If more complex curves 
are needed to describe the structure, or if there are 
Eusocontinuities, additional SEL's сап be used to meet the 


requirements of the particular use. 
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Coordinate transfer with 


Nodes repositioned equal increments 
exponentially . 





Coordinate transfer with 
unequal increments 


Fig. 10. Mesh generation with and without equal element size 
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TABLE IV 


ISOPARAMETRIC ELEMENT TYPES USED 



















Common 
Name 


Number of 
Nodes 


Exactly 
Represents 





Linear Four 


Straight Line 





Quadratic Eight Parabola 





Cubic Twelve CUBIC  ROlIynomia 


Шш ОООО п Go note that SEL boundaries which 
cannot be exactly described by the element type chosen (i.e., 
are of higher order than the element used to describe them) 
Br eome sliehtly distorted in the mapping procedure. Except 
Клозе boundaries Which coincide with the boundaries of 
EN trüctbture, this is of no consequence. The structure 
boundary nodes, however, must be corrected to locate them 
properly in the real coordinate system. A provision to 
accomplish this for the curved boundary of Face B (Fig. 9) 
has been included in PLIMEG. 

A representative mesh, as used in this study, is 
shown in Fig. 11. Quadratic elements are used througout 


mes Study. 
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SE ondary Condition Calculations 

I mec Ncondgltvelons tor the structure are in two 
categories — force and displacement. As established in 
section III D , there are two sets of force conditions and 
αν Sets Of deflection conditions. These are easily coded 
into PLIMEG where the appropriate input data for the FEM 
program is prepared. A particular set of boundary conditions 
is imposed by a code number which exactly corresponds to one 


ШОШ е four problems Outlined in Section III B. 


D. SUBROUTINE PROBLM 
mime Gesu Structure, divided into the SEL arrangement 
ШЕ ed. is illustrated in Fig. 12. SEL numbers and SEL node 
numbers are indicated. The SEL arrangement was selected 
after considerable experimentation with other arrangements, 
based on the following considerations: 
l. Small elements are required at the base of the notch 


because of the high stress gradient there. 


2. Large elements are adequate away from the notch, 
specifically at the end of the structure where loads 
are applied. 


9S. Storage capacity of the computer limits maximum number 
of elements to not more than 180. 


6. Extreme element corner angles (i.e. near 0° or 180°) 
Must be avoided. 


pe Each SEL must be a four-sided figure. 

6. Relying on the principle of Saint-Venant, a single 
element boundary may be used to represent the loaded 
end. 

The coordinates of most of the nodes are obvious after 


certain key nodes are located. The origin is at node 27. 


Nodes 1, 4, 7, 20, and 33 are established by the parameters 
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directly (R/W, D/W, 6), when L = 5W is provided. The x 
eoordinate of node 30 was established at W+R and the y 
coordinate of node 14 at D/6. Boundary A is an are of the 
parabola passing through nodes 14 and 20 with zero slope at 
node 14. Boundary B is an arce of another parabola through 
nodes 4 and 17 perpendicular to boundary A at node 17. 
Boundaries C and D are straight lines. Intermediate nodes, 
where indicated, are located at the 1/3, 2/3 points along 
the curve establishing the boundary. Table V is a summary 
the SEL node coordinates. 

Foesprocrem PROBLM, written to calculate SEL node 
coordinates uses the above criteria and requires a single 


card input: parameters R/W, D/W and ө. 


КОО SUBROUTINE ANS 

Mae answer required for each problem is the stress 
Bemcentration factor which is calculated .in subroutine ANS. 
This program produces a single card reporting problem 
Parameters and the stress concentration factor associated 
with them. The stress concentration factors reported are 
based on the following definitions of ο σαν 


fase 1, uniform loading: 


“hom s Eu 


where: 


I. 


total load applied * 
A = 2D = net width of the bar 
*Actually here and elsewnere, loads are per unit thicxness. 
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TABLE V 


SEL Node Co TIT ces 


y coordinate 













ELS x coordinate 
d 0 D 
1 Ευ 5) Note (1) 
1 К sin(60° - 28/3) Note (1) 
NS R cos 6 . (DER) - R sin Ө | 
J xy + (W-y)) cot 0 W 
5 Xp + (7 77Xp 2/3 Note (3) 
7 xy  2(X 7*23 Note (3) 
lo 0 D/6 
1,2,5,6 R+W Note (2) 
50 хут і (x497X4 42/3 Note (2) 
5,6 X17 * 2(x5974 42/3 Note (2) | 
56 DW W 
2 0 0 
2,6 o | 
| 
IE 


Notes: 


EB NDetermined by circular arc of radius R through 1,1. 


(2) Determined by parabola A through 14, 20 with zero slope 
at 14. 


(3) Determined by parabola B through 4, 17 perpendicular to 
parabola A at 17. 


MPoints not listed are on straight line boundaries equally 
spaced between corner nodes. 
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Case 2, bending moment: 


g - Mc/I 
nom 


where: 


M = total moment applied 


ο DP- distance from neutral axis 


I 2D3/3 = moment of inertia 


esse 3, uniform loading: 


nom = P/A 
where: 
P = total load applied 
A = D = net width of the bar 


Case 4, bending moment 


σ = 

vn Me /I 

where: 

M = total moment applied 

ce = D/2 = distance from neutral axis 
I = p2/12 = moment ofl inertia 


F, PROGRAM MAIN 
The four programs: PROBLM (establishes SEL node coor- 
dinates based on parameters), PLIMEG (discretises the struc- 


ture and computes boundary conditons), PLISOP (solves the 


! 





stress problem) and ANS (calculates the stress concentracion 
factor) are joined together by MAIN. All card output, except 
une single card from subroutine ANS, is suppressed. Input 
data and intermediate results are passed within the program 
common blocks and via computer direct-access devices. 

The functions of the main program are to read and store the 
mae data, call the various subroutines in turn, provide 
ОШООО of the solution cards and prints, and to increment 

the problem parameters so as to run through a related series 


er evaluations. 





V. PROGRAM OPERATION 


Е INPUT 

The problem input data deck consists of nine cards, a 
Mic al example of which is shown in Fig. 13. Card 1 lists 
Pe problem parameters: R/W, D/W, and 6. Cards 2 and 3 
are title cards which may contain any message to identify 
the run. Card 4 lists the mesh parameters: number of SEL's, 
number of nodes in each element (not SEL), problem τοπ. 
eon number and an indicator for plotting the mesh. Cards 
5 through 8, are SEL connectivity cards. Each card lists 
Bae number of increments to be used in the ë and n directions, 
reir” respective exponential coefficients, and a listing of 
the nodes, appropriately ordered, defining ne corners cf 
ο ου Card 9 lists the solution parameters: number of 
Esas, an indicator for the plane stress solution, and 
ode or the numerical integration technique. 

The job control language cards (JCL) needed to make the 
C; Operational, details of card format, and program 


F Ons are discussed in Appendix D. 


Ee OUTPUT 

micnough the output required in this particular applica- 
tion is a single number (SCF), a οι more is available. 
All of the input data is echo checked and the following 


important data passed between subroutines is printed: 
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13 SAMPLE INPUT DECK 


FIG. 





ὑπ ποσα Coordinates. 


Element node coordinates with the load boundary conditions 
as imposed. 


Element connectivity. 

Displacement boundary conditions. 
Material properties. | 
Displacements u and v for each node. 


Reactions at the nodes where displacement boundary 
conditions were imposed. 


A Torce equilibrium check which is the result of summing 
all x-direction and all y-direction reactions. 


Average stress components at each node. 
Average strain components at each node. 


otress concentration factor. 


Although the stress concentration factor is independent 
ο. oung's modulus and Poisson's ratio, the values 30 x 10° 
and 0.3 respectively were coded into the program as represen- 
Eve of common engineering materials. The average stress 
and strain components are computed and printed in two ways: 
the cartesian coordinate components rohe and e) and the 
principal values (04,05, and or) гони спе inclination orf 
the principal axes (4$). 


ener output in the form of data cards and the plotted 


mesh has already been mentioned. 
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Vive: ROGRAM OUTPUT ANALYSIS 


A. EQUILIBRIUM 

static equilibrium is maintained in the finite element 
Lion of all e studied. As indicated earlier, the 
equilibrium of forces is quickly checked from program output 
by comparing the imposed loading with the reaction summations 
printed by PLISOP as "equilibrium check". 

mime Check On moment Summation is slightly less convenient. 
meeeran output, including loads and reactions at all nodes, 
апа the coordinates of nodes, is available but the actual 
arithmetic must be done by hand. These external calculations 
were made for several representative problems. The results 


Were as required: moments and forces sum to zero. 


B. MAXIMUM STRESS LOCATION 

sane tre program output, the location of the maximum 
principal stress was determined for several parameter combina- 
tions. In every case the maximum occurred tangent to the 


meee at the bottom of the notch. 


EE CONVERGENCE 

To investigate convergence behavior, we first identify 
το S and columns of the SEL divisions used. The columns 
EXE referred to as band A and band B; the rows as band с апа 
band D (Fig. 15). In what follows, band D always consists 


mone row of elements. This limit was imposed because 
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Br core storage availability in the computer used. 

The refinements were accomplished by first assigning a base 
Ember to each of the other bands, then multiplying by a 
Successively larger integer for each new mesh. Table VI 
shows the number of elements per band and total elements in 


the mesh. The mesh ultimately used is included for comparison. 


TABLE VI 


Convergence Subdivisions 


Ne otal = (N fN) (N +N.) Meter co mie 14) 


Band First Second Phra Fourth 
Mvision Run Run Run Run 
(1) (2) (2) (29, 


D 


ο στα] 
Elements 





(1) This division was used primarily to verify correct 
performance of program. 

(2) These divisions were used to verify aná evaluate convergence. 

(3) This division was used for volume data generation as an 
optimal compromise giving good accuracy and Satifactory 
computer turnaround. 


D 





The refinement was dorie for several combinations of the 
problem parameters. Figure 15 is a representative plot of 
the stress concentration factor vs. (reciprocal of) total 
elements in the mesh for two combinations of the parameters 
as shown. In both cases only the three finest meshes used 
in the convergence study (50, 105 and 180 elements)are 
η. Fitting a straight line through the points and 
wapo lating back tothe axis of ordimrates, an approximation 
for the converged result can be obtained. The value obtained 
using 135 elements (see below) has been added. Note that it 
plots on this same straight line. 

Consider the upper curve shown which represents the case 
of a shallow notch. The extrapolated value for stress 
concentration factor is SCF = 1.773. The 180 element result 
is less than 0.8% high while the 135 element mesh used is 
less than 1% high. Results from studies using other parameter 
combinations, including 8 # 0, showed the same convergence 
pattern. 

The mesh of 135 elements was chosen for data runs because 
that size mesh requires approximately 3/4 of the computer core 
Storage space available. From a practical standpoint data 
collection proceeded without undue special arrangements in 
the computer center. The larger number of elements required 
elaborate special consideration by operating personnel which 


was not warranted by the 0.2% difference in the results. 
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Fig. 15. Convergence Plot 





D. CORRELATIONS 

Nerds ta collected dls Cisecussed in the next section. As 
will be seen, experimental data from the literature correlates 
very well with values of SCF calculated. The trends established 
are verified and the experimental values are sometimes 


identical to three significant figures. 
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ο DIC CUSSION OF RESULTS 


A. RESULTS PRESENTED 

The SCF's presented in Appendix A document a represen- 
Drive variation of the parameters for problem 1. It is 
EBsvsous that either 0 - 90? or D/W = 1 constitutes a bar 
Ea mo noteh and results in Stress concentration factor = l. 
Several computer runs with D/W very small showed that as D/W 
approaches zero, SCF indeed did approach unity. 

To establish the curves, D/W = 0.1 was taken as the 
initial value and incremented by adding 0.2 through 
D/W = 0.9. Likewise, R/W was incremented by steps of 0.05 
πο. through 0.5. 8 was incremented by steps of 15? 
com 0% through 60°. 

The data are presented as several families of curves. 
Each family has a specified 8 associated with it. The curves 
are plotted as SCF vs. D/W with R/W as a parameter for each 
family. Linear interpolation between the curves may be used 
Without appreciable error. Interpolation between the families 
for 8 = 60% and 0 = 90%, where stress concentration factor 
one, Should not be attempted. Data points actually 
calculated are indicated by a cross. The curves through 
16 points were generated by fitting a sixth order poly- 
memial through them. 

Αν ος тїп was Made to prove the programing for problem 
2; no data for this case are presented. One run was also made 


for each of problem 3 and 4. The data collected are presented 





ШОШ ес оп VII Din the form of a correlation with experimental 
usta from the literature. 

Time was not available for more extensive data collection. 
However, now that the software has been perfected and proved, 
NONIS a perfectly straightforward task to obtain additional 
data as computer time becomes available. Appendix D suggests 


the parameter ranges which should be investigated. 


E EFFECT OF NOTCH DEPTH 

Selected data from the studies of problem 1 are plotted 
in Fig. 16. The parameter D/R is constant. Experimental 
data plotted is from Flynn and Roll; however Kikukawa's data 
shows the same trends. The data in the experiments was for 
the symmetrically located U-shaped notches as in problem 1, | 
with 6 = 0. .Both Neuber's results and Flynn and Roll's 
experimental data are included for the hyperbolic shaped 
symmetric notches with the same parameters. 

The general close agreement between the FEM and experimental 
values for the U-shaped notch is easily seen. The largest 
Gifference noted is for W/D = 4.33 where the difference is 
less than 2%. 

eé micht examine the difference by utilizing an aspect 
ratio defined as effective length from the notch to point of 
‘load Es cPemnsgsvded by notch depth. Table Vel compares 
Meese Value. AS can be seen, the calculations alisw a much 
more favorable (i.e., larger, and thus more realistic) ratio 


pan that provided in the experiment. 








1 2 3 4 W/D 


Fig. 16. Effect of notch depth on SCF for RAW = 1.5 


Q——————O Numerical U-notch (Finite Element) 
Ε---------- Experimental U-notch (Flynn & Roll) 

$—— — — —Q Experimental Hyperbolic notch (Flynn & Roll) 
————— — Analytical Hyperbolic notch (Neuber) 
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TABLE VII 


Aspect Ratio L/(W-D) 


FEM 10 fees. ο. 
Calculations 


Experiments 8.6 13 2 wO 





The maximum point on the curve in Fig. 16 occurs at 
W/D = 2, a relatively shallow notch. Experimental data 
prehes at this point exactly. Note, the error that would be 
made if Neuber's hyperbolic notch result was used for such a 
notch is more than 7%. Experimental data presented by Flynn 


and Roll point to a systematic error in Neuber's results. 


C. EFFECT OF ANGLE 8 

The data plotted in Fig. 17 compares the stress concentra- 
fon factor for 8 = 60° (SCF eq) to the stress concentration 
factor for 6 = 0° (5080). The parameter used is W/D. Two 
experimentally determined curves show that increasing the 
angle of the V from 0° to 60° decreases the SCF cq. This 
effect is more profound as the WD ratio is increased (i.e.; 
for deeper notches). 

ΙΙ; chescata as in Fig. 16 shows the variation of 


SCH 0 with changing W/D for constant SCF Note the close 


0° 


wne lation of the experimental data and FEM data. 
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INL. Effect of V angle on SCF. The dot-dash line 
(SCF, = 5080) is included for reference to emphasize that 


SCH zo is always less than SCH, for any value of the other 


Merameters. Similar data could be plotted for other values 
oe О, 
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E W/D 
foe. 18. Effect of W/D on ΒΟ ΚΩ these lines 
represent profiles (for SCF, = 2.3 and SCF) = ο ο. οἳ “the 


data shown in Fig. 17. The same plot symbols are used as 
in Fig. 17. | 
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Figure 19 demonstrates that increasing the angle Ө 
meouces the stress intensification and this effect is more 
Profound for smaller radius notches than for notches where 
the bottom radius is large. That is, the SCF for a given 
notch may be reduced by increasing either R or 0 while the 


other parameters are unchanged. 





W/R 


πο. Effect of angle 6 on SCF as function of radius 
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Dr SINGLE NOTCH GEOMETRY 

The only experimental data available dealing with the 
single notch geometry is that presented by Cole and Brown [5]. 
ST data is for a particular ease of loading such that the 
Meme Of action of the applied force is through the midpoint 
ШИ еспе width of the bar. By an appropriate combination of 
the SCF's calculated for problems 3 and 4 this loading can be 
reproduced as shown in Fig. 20. The defining relationship 


for SCF as used by Cole and Brown is: 


σ 
_ "max 
ΠΠ 
where 
ux - maximum stress 
P = total load applied 
D = net width 


The defining relationships for problem 3 (SCF,) and problem 
Ц (SCF, ) are contained in Section IV E. 

For the parameter values D/W = 0.5, R/W = 0.1 and 6 = 0°, 
Brand Brown report SCH, = 1.95. The computed results are 
SOF, = 7.85, SCF, = 1.08. Makine the appropriate substitutions 


Solving for maximum stress: 


62 








= 


το υπυ load (P) 





Uniform load (o = ΡΛ) 


Bending moment [M = (W-D)P/2] 


το. ο. Bar with one notch ον Ppoint oad decomposed 
into constant and linearly varying components. St. Venant's 


principle accounts for the absence of higher order 
components. 





Q 
| 


P ‚Mc 
max - SCF 3(5) SCF) (=>) 


| P 
(SCF, η SCF} ) (5) 


1.91(Ë) SCF(5) 

Thus the FEM analysis gives SCF = 1.91 compared to Cole and 
Brown's 1.95. There is approximately 2% difference in the 
results. A special computer run with this specific loading 
Ane SCF = 1.91 also. This calculation Served to verify 
meer pOsability of the two cases and the correctness and 


Euvernal consistency of the programming and the inputs. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The FEM program as presented in this thesis can be used 
posesteulate the stress concentration factor for the specific 
problems listed, over a wide range of the parameters used. 
Further, by appropriate superposing of these four specific 
Pmootems, the oUP can be calculated for bars which are notched 
either on one side only or symmetrically on both sides and 
are subjected to general tensile and bending loads. 

One Of the conclusions from an analysis of the data 
obtained for problem 1, the case of the symmetrically notched 
Ber subjected to a tensile load only, and a comparison with 
pabulations and recommendations in the literature is that 
the three parameters (minimum radius of curvature, gross 
width of the bar, and net width of the bar) do a rather good 
job of characterizing the geometry of the notch, in the sense 
that interpolated values of the SCF obtained from available 
data agree rather well with the more accurate values obtained 
iy FEM analysis in this study. 

However, the analysis also shows that there are some 
small but systematic deviations that are assignable to the 
fact that these three parameters are not sufficient to 
characterize the geometry of the notch completely. The data 
puesented in Appendix A accordingly provides a more complete 
ENUusceurete representation of SCF's in this case than was 


previously available. 





The stress concentration factors presented in the litera- 
ture for several specific cases of V=eMapedeand Vasmaped 
notches are essentially correct. However, the data presented 
NS Iimited, mainly because of the difficulties associated 
with experimental investigations. This program overcomes 
those difficulties and can solve for a single SCF in less 
than three minutes of computer time as compared to several 


Nours by accepted "experimental" methods. 


B. RECOMMENDATIONS 

Becordineily,s at Ts recommended that the program described 
herein be used to determine SCF's for problems 1 through 4 
and that the results be published for general use. Appendix 
D presents a specific plan for obtaining the necessary results. 

Buecher, IT interest in other notch shapes warrants it, 
MS program should be appropriately modified and the SCF's 
puoeulated and published. 

An axi-symmetric stress analysis program has been written 
SES Eahkshandehpour [3| as part of his NPS thesis. That 
meorram could be similarly adapted to calculate SCF's for 


general axisymmetric geometries having notches. 
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APPENDIX A 


CURVES OF RESULTS 


The stress concentrations presented in this appendix are 
for the problem of symmetrically located, round bottom, 
V-shaped notches subjected to uniform tensile loading. The 
Gata collected is plotted “as crosses. A sixth order 
polynomial is shown fitted through these points. 

' The curves are plotted as SCF vs. D/W using R/W as a 


parameter with each family representing a single 6 value, 


where 
lau i οἳ the bottom οἱ the notch 
D = net half width of the bar 
W = gross half width of the bar 
6 = half included angle of the V 
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Stress Concentration Factor (SCF) 


p indicates FEM data. 


σ 


max 
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APPENDIX B 


NEUBER'S NOTCH ANALYSIS 


This appendix is included since we have been unable to 
find a reference which quotes these useful and important 
Mmeewlts in à brief and clear fashion. 

We are concerned with the two-dimensional case of an 
axial load transmitted in the y-direction through the net 
section of a thin bar (plane stress case) having a hyperbolic 
ШО ету. All quantities which follow are for unit thickness. 

To describe the geometry and to provide a basis for the 
ού του for elastic stress distribution, Neuber employed 


the mapping 


c cosh E, Е = cosh" t (z/c) - in[(z/c)*(z?/c? - 1)$j 


2, = 
(Bla,b) 

meere z and ζ are complex variables 

ο στ κι iy , ¿£ = £ + in τα ο) 
ele i is the imaginary unit. In real terms this results in 

x= cC cosh — cos тп; у = c sinh E sin n (Bsa. by 

E 1 -1 7ك‎ 7 -1 

E = %cosh "(GtH) , n= cos [G-H] (Blla,b) 

where 
> 
ee ae « 2(y9-x^)/e? l (B5a,b) 
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This mapping provides a correspondence between the xy 
cartesian system of coordinates and a &n system useful for 
eining the geometry of the notch. Lines of constant E are 
ellipses 


( X y? 


as 
c cosh Ẹ E τ 


B E: ΠΠ ἐ 


tte Nines Of constant n are the hyperbdolas 


c ^ 1 (BT) 


Figure B-1 shows the pattern of the En network. 


š 


E = constant 


= constant 
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mee. BL. Neubder's hyperbolic, symmetric notches in an 
infinite plate. 
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A complex Airy stress function is used to solve the 


EN ουἱ1ο τν problem so that 
сыз т = 0 (В8) 


along the notch boundary which is the hyperbola 


m y^ cot np = c° cos “no | (Β9) 


eommesponding to the value n = no: The half width at the 
base of the notch is the value of x when y = 0. We denote 


this by a, so that we are able to evaluate 
ο ο εστι. کک‎ (B10) 


This choice of c selects a definite shape, namely that 
which has net width 2a and whose hyperbolic boundary has 


asymptotes 
ο Cate (B11) 


IAS significant to note that except for a change of overall 
scale, the value of No is the single parameter which defines 
the geometry of the notch. The two "halves" of the hyperbolic 
boundary are fixed relative to each other by the fact that 
their mutual asymptotes intersect at [πη opigin x = y = 0. 

One of the most significant derived quantities is the 
Eon radius of curvature of the hyperbola, p, which occurs 


at the base of the notch. 
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p=a tann, : 2) 


RE oO on που given in detail here, permits finding 


the Stress state (0,,0 0 xy) or E n at al, Point orf 


y 
ο. τιοῖσα. A constant, 0, is involved in this solution. For 


our purposes we need only the stresses along the line 


T= O, Viz: 
= = 3 Е 
с σῃ ¿scan sinn sin no) 
= O, = 2b 3 (si en + sin“ ) Came be) 
ο ο. Csen sinn sinn, 250 De 
= = 0 
Тху "En 


One of the results we will use in Appendix C is the sum 
DU UE lb ese n | (B14) 
along y = 0. 
To determine the constant b, we note that the tensile 
force (per unit thickness) transmitted in the y-direction is 
EE T2 J б. ах (5. 
0 


De titutins Eq. (Bl3b) and solving for b, we get 


О КГР cos Ng)/2a(m - 2ng + sin 2ng) (B16) 
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ИО Е Ба (B16) into (B13b), and evaluating at 
n * ng, we get the maximum tensile stress (at the base of 


the notch) to be 


ON E (2P cot по) /а(т - 2ng + sin 2ng) (B17) 


The nominal tensile stress is 


P 
Onom 2a (B18) 
Eud the stress concentration factor is 
SCF = 0 /σ = (4 cot Ng)/(m - 2ng + sin 2ng) (БКО) 


max nom 


In terms of the ratio of minimum radius of curvature p and 


net half-width a, we have 


men = 4/kin = 2 arctan k + 2k/(1+k°)] (B20) 
where 


К = р/а (B21) 


Mhis shows the SCF in terms of the ratio of two evidently 


significant and measurable geometric dimensions. ' 
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APPENDIX C 


THREE DIMENSIONAL EFFECT 


ШОО PREAMBLE 

merore proceeding with this discussion, the following 
points should be made clear. First, there is every reason 
to place full reliance upon the accuracy of the determina- 
tions made by the two-dimensional FEM program, that is, to 
assert that these calculations accurately represent the 
theoretical two dimensional problem (plane stress, as 
presented here). From this standpoint, the study of three 
Gimensionality is intended only to assess the valididty of 


photoelastic determinations which were compared to the FEM 


determinations, and we in no way depend upon three dimensional 


Bem@avior in presenting, with great confidence, our results 
for the two dimensional case. However, one cannot evade 


the responsibility of making some sort of evaluation 


of the effect of finite thickness upon determinations valid, 


ШОО УОС Ту Speaking, only for cases of zero thickness or 


Иа се thickenss. 
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DA SPECIFIC PROBLEM 

In addition to other difficulties of photoelastic evalu- 
... τ. there is an inherent error which is due to the fact 
that all photoelastic models have finite thickness so that 
the conditions of plane stress are violated, at least to some 
rent. 

MOS first approximation, as shown by Timoshenko and 
Geodier [15], the variation of a опаа stress field, 
due to finite thickness, may be obtained by employing, 


instead of the usual Airy stress function Py (xy), satisfying 











y 3^ 9 2 0 И? 
V M CMS a σν 25: шо = o. Я = Ту (Ca p O n) 
dx ΟΥ ox oy 
σαι τε function 
S 
$(x,y;2) 7 $9 G3) - 503; 9900) (C2) 
where 2 2 
dx ΟΥ 
Mie (modified) stresses 
3* 32 e 
ج‎ RA : ο. — s DM ok (CHa BIC) 
ΟΥ dX дхду 


mead irom this modified stress furction, thus have a 


MemebOlic dependence upon the third dimension, and, in general, 


19 





agree with the elementary (strictly two-dimensional) solution 
Epub on the plane of symmetry, z= 0. 


νο ους TE the first tensor invariant of the 


0 


elementary solution, that is 


“о (σ.})ᾳ = ος (C5) 


where we have added the subscript 0 to indicate that these are 
the elementary or unmodified stresses. 

In photoelastic experiments, the light beam is affected by 
phe distance-average stress-state and this differs somewhat 
luu the stress-state at mid thickness. It is the latter, 
however, which corresponds to the theoretical plane stress 
distribution which the photoelastic evaluation is intended 
Borreveal. 

determine tbe correction it is necessary to have, 

ШЕЕ Су. second derivatives of the first tensor invariant. 
Meese involves two-fold differentiation of a stress field 
which itself has been determined, usually, by graphical or 
numerical integration procedures. The problem is somewhat 
easier to deal with at points on the free edge. Here there 
v One stress, namely that in a direction parallel to 
de free edge, with which we are concerned, and all that is 
required is two-fold differentiation of 9 inta direction 
normal to the edge. Nevertheless, it is quite unlikely that 
sufficiently accurate stress determinations will have been 
Mame Orr the interior of the model area to permit а reliable 


evaluation of the second derivative needed. 
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In the particular case with which we are concerned in 
ΙΙ thesis, namely the matter of stress concentrations at 
the bottom of notches, we can construct a plausible theoretical 
means of obtaining the desired information. This procedure, 
and experimentally accomplished equivalent procedures are, 
however, of doubtful accuracy and probably do no more than 
glve a very rough estimate of the proper correction which, 


Memyuriately, usually is not large. 
2 


jl problem, specifically, is to determine : 99 
dn* 
where 99 rove oUlimOm principal stresses πα αι πο σας 


direction normal to the boundary at the point at which the 


correction is to be made. The evaluation is at the boundary. 


BRZZEROCEDURES TO SOLVE THE PROBLEM 

The following procedures suggest themselves. First, one 
ШЕШ πε the photoelastic stress information, properly integrated 
so as to obtain the stress-sum, and perform two-fold 
differentiation. Or, still proceeding experimentally, one can 
norm another experiment, using an interferometric procedure, 
to get the isopachic lines, which provide a direct experimental 
evaluation of the (average) stress-sum, and then differentiate 
m еа twice in the normal direction. Third, one can 
analytically determine the desired second derivatives by 
moving ah available theoretically exact solution to a 
geometry which "matches" that of the case under investigation. 
ΠΕ, One can construct a photoelastic model of such a 


geometry for which an analytic solution is available and infer 


eal 





the desired correction by comparing photoelastic results 
with the theoretical results and then apply the same correction 
to the case under investigation. 

Although superficially these procedures seem reasonable, 
each has flaws of fatal or near-fatal magnitude. The first 
two methods described do indeed relate to the geometry of 
Mmmemease at hand rather than to some substitute geometry. 
ШОС in theory, they are capable of providing the desired 
mun Lion. However, the fact that a stress-sum field must 
be constructed with sufficient accuracy that a second 
derivative can be relied upon indicates the need for a degree 
of precision considerably beyond what might be reparded as 
normally achievable. 

The third and fourth methods are essentially the same — 
one is dealing with a substitute geometry for which a 
Bneoretica] result is available, In the third method one 
needs a more detailed mathematical representation of the 
Scion; in the fourth method all that is required is the 
theoretical value of the edge stress — and a second photo- 
elastic determination. In each of these methods however, 
the reliance is upon the assumption that the second derivative 
mu he stress-sum is the same for the substitute geometry as 
u Case of interest, or at least that the ratio of this 
second derivative to the edge stress itself is the same for 
Ene TWO cases. 

τ. τ: quite disinelined to believe that there is any real 


Mita rFicatíion for such reliance. If the substitute geometry 


82 





(with known analytic solution) is sufficiently close to the 
actual geometry that one may rely upon the coincidence of 

the second derivative of the stress sum, then indeed there 
Should be even closer coincidence between the edge stresses 
Fu C ves so that a photoelastic investigation is not called 
for at all. It is the very fact that the geometry of interest 
does not coincide with a geometry for which there is an 
available analytic solution that causes us to make the 
mootelastic study in the first place. 

However, even though we may have deep and quite reasonable 
bts concerning the accuracy of our evaluation or estimate 
of the second derivative of the stress-sum (or its equivalence 
in other terms), nevertheless, for КОО КО О RET models 
ОШОГО is reason to believe that the corresponding variation 
in edge stress through the thickness of the model is "small" 
so that a greater degree of error may be tolerated in this 
estimate than would be tolerable in the basic evaluation of 
the edge stress. 

Our own study of this matter has involved (1) analytical 
differentiation of the stress-sum field for a known geometry 
Mabproximating" that of the case under study, (2) numerical 
differentiation of the stress-sum field cbtained by use of 
a two-dimensional FEM program, and (3) actual use of a three- 


dimensional FEM program. 
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In the following paragraphs we concern ourselves with a 
U-shaped notch. The radius of the bottom of the notch is 
1.5 inch, the gross half width of the bar is W = 1.3 
Baches, the net half width is D = 0.3 inch and the half 
Inuekness 15 0.125 inch. The angle 0, as used in the text, 
EO = 0°, 

a. Analytical Differentiation 

1. Neuber's Deep Notch 
The geometry of Neuber's hyperbolic notch in 
Eununisnite piste is shown in Fig. Bl. Equations B13 and 
B14 in Appendix B give the following expressions for 


stress: 


Ge) 5 csc?n (sin^n * sin?n.) (СО) 
yO O 
66 = (σ * 0. )p - lb ese n (C7) 
Letting csc N = (1 -u2)717¢ апа designating t. DV ^ ux 


we get 


-5/2 ese“n,/D“ (C8) 


0" = дь (142u°)(1-u°) 


Perens y = 0, 0. varies both with respect to x and z; letting 


c = v/2(1+v) ; 


04 





we get from eq. (C4): 


Oy = (σι) ο) - ο 6" 2 ССО) 


Maine 2A = t , the thickness, and evaluating the constant, 


D: 
р А 
p = f f us dz dx (C11) 
(0 Ὁ 
b = p eos No/HADEG - 3H/2] (C12) 
where G = n-2n,*sin 2ng (C13) 
Ἡ - λε ц у 
= Ас cot τ, (CI) 
Eom eq. (C10) 
ο... = Hb ese No (ORO 
2 
E Е 11 
er η c 6g" A (016) 
Ше сале 
Ао, = Care - КЕ (CP 
Ac. 8) 
m ES H ese n, (142 cos na) τ 
max 0 








Mem Cle πες Profile through the thickness 


Figure Cl shows the stress profile through the model thickness. 
The photoelastician sees ©., averaged through the thickness 


y 
2А: 


Ас 
) = (0) Е ск (с19) 


σ σ 
ож ауе y max 


( 


wider the case: D = 0.3, p= 0.2, A = 0.125, 


p = D tanón, . We find 
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ο 278285 


zn = 0,982 O P naz ο 205 


( 


(Here and later 90. GenoOtescuvie evaluation of 89" ac. the 


base of the notch.) 


Maat is to say the photoelastic evaluation of stress provides 
an uncorrected value for ee that is nearly 2% low. 
2. Circular Hole in a Plate 
One substitute geometry which suggests itself 
ås that of a circular hole in a plate under uniaxial tension. 


Figure C2 shows a round hole in a large plate. 





ОШКО МЕ тоссо е in an infinite plate 
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If the plate is large in comparison with the size of a 


peca lar поте the following relations hold: 








o = S[(1-R®) + (143R"-4R?) cos 20/2 (C21) 
E 2 lj 
σα = S[(1+R ) - (1t3R )ocos 2a]/?2 (C22) 
θρ pes т a Sl ARCOS ρα) 023) 
where R = ro/r (02:5 
Introducing the three-dimensional behavior, we 
have 
2 2 
2 340 340 
(с) = 52 = — `Š + ezf گی‎ (C25) 
or or ar 


and evaluating at the edge (r = Po?» we get 


E ace M ος 


[S-2S cos 2a] + [6S (z/rg)* v cos 2a/(1+v)] (C26) 
iabilterclis nos tevatuations of this formula 


Оо points around the hole, in terms of the ratio 


λ/τρ where A is the half thickness. 
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TABLE C1 


EDGE STRESS AROUND HOLE. 





Inserting the numerical data as previously used: r, = 0.2, 


feo. 1253; and solving at a = 90°: 


(с) 


edge 50 ση (C27) 


Proceeding as before: 


(σ) = 0.91 (о) (C28) 
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A more precise analysis using Howland's [8] 
ulation yields essentially the same result. It is apparent 
It a circular hole is used to infer the needed correction, 
as suggested in Section 3 above, SCF's reported for the 
U-shaped notch will be slightly high. This is apparently the 


procedure used by Flynn and Roll in their experiments. 


4. NUMERICAL DIFFERENTIATION FROM 2D FEM 

For this evaluation, we deal directly with the U-shaped 
morcn, The same numerical data was used to facilitate direct 
comparisons. Our objective is to take the FEM computed values 
of Dx + Sy Io cesa tone the line normal to the notch, 
fit a polynomial through the values obtained and compute the 
End derivative in the normal direction. Figure C3 shows 


Bue notch. 








EE ς 


Fig. C3. Symmetric U-shaped notches (not to scale) 


90 





Examination of the stress-sum profile reveals that a 
moot Curve is not obtained even with the finest mesh used. 
(This is primarily because of a minor systematic inherent 
Giscrepancy between the FEM evaluation of stress at a midside 
node and at a corner node.) Choosing the element corner nodes 
Esa fourth degree, least square polynomial ean be calculated. 
Examination of the second derivative of this polynomial reveals 
fiero tt TS sTowYy ”“converftne"to"”some “value with each 
Successive refinement of the mesh. Viewing each of these 
ws as a term jn a convergent sequence, a convergence- 
accelerating procedure due to Shanks [14] was used to estimate 


09" ο ὃς base of the notch. The numerical result was: 


(Gis) ο ο. 5.) (C29) 


X ave x max 


bee DIRECT COMPARISON FROM 3D FEM 

The FEM program used in this study was TRISOP [10]. 
Quadratic elements were used for the solution. The mesh was 
generated by TRIMEG [1], modified in the same way as the two 
dimension mesh generator (PLIMEG), to progressively vary the 
element size. The mesh in the x-y plane was the same as that 
r (in the 2D solutions. Two equal size elements were used 
he half thickness in the z direction. Because of symmetry 
Beemer che midthickness plane, it was sufficient to consider 
ami half the plate, The half thickness of the bar (previously 


designated A) was 0.125. Figure C! shows the notch. 


2 








Fig. C4. U-shaped notches. (xy, yz, zx planes are planes 
of symmetry. The planes z = +A are outside faces, 
the planes y = #w are the unnotched boundaries) 
Each computational run involved nearly three and one-half 
hours of computer time. Because of on a few runs were 
made to insure proper program operation and gather this date. 
Since a large data gathering task was not contemplated here, 
ШО sh generator and FEM program were not coupled together. 
All data input to the mesh generator was manually prepared. 
All calculations subsequent to determination of nodal stresses 
Was done outside the computer. 
A parabola was fitted through the values oe at the corner 


Berges along the line y = a = 0.3. The result of this Procedi p 
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cast again in familiar terms was 


$9" = 14.98 


(o_) = 0.991 (o,) (C30) 


x ave X max 


The former value was inferred from the coefficients in the 


parabolic fit. 


oe) CONCLUSION 

We are now in a position to compare the various results 
and estimate the importance of the three-dimensional effect. 
The aspect ratio of the bar Тоос саа һеге 15 


thickness ` 0.25 . 


Dueh ome ο τα (C31) 


The results obtained are Summarized in Table C2. 


TABLE C2 


Results of Three-Dimensional Analyses 


Method BO Remarks 


Neuber notch ο 2 18.85 Two derivatives of 
analytical stress-sun. 


Hole in plate 0.94 99.79 Two derivatives of 
analytical stress-sun. 


2D FEM 0.994 10 σος τι αρ, τος οὗ Έῃδ 
curve Lit through 
stress-sum from 
numerical data. 


3D FEM 0.991 14.98 Parapolası itstchroush 





numerical computations 
CRO 
x 
== we] 








As can be Seen from the results presented here, the 
effect of three-dimensionality is indeed small. The two- 
dimensional and three-dimensional FEM studies compare 
closely enough to permit extracting the third dimension 
effect from PLISOP. The extraction is not easy. However, 
the procedure outlined is readily accomplished with the 
computer and requires much less computational time than does 
mere exectutaron “of DRISOP, 

Mersee also irem these results that an indiscriminate 
adjustment of photoelastically obtained values of SCF based 
on the SCF for a known geometry may cause overcorrection. 
The overcorrection is small, however, and the error causes 
conservative estimates, that is, one overestimates the 


rL JO (o σ Р 
- ( d cU Boh 
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APPENDIX D 


SUGGESTIONS FOR PRODUCING ADDITIONAL NUMERICAL RESULTS 


press cOnecenyration factors for the four problems 
E Led in this thesis can easily be calculated by using 
the program in Appendix F. For all four problems the useful 
range of the parameter D/W is from zero to one and the range 
of 0 is from 0° to 90°. Note that either D/W = 1 or 9 = 905 
КО bes a bar with no notch and results in SCF = 1. 
D/W may be incremented in steps of 0.2 (program statements 
110 and AA 01200) and 8 in steps of 15° (program 
Statements AA 01380 and AA01180) for all problems to provide 
sufficient data. The incrementing step selected in the 
program for R/W is 0.05 (program statements AA 01360 and 
AA 01190). 

The SCF's for problem 1 over the full range of D/W, all 
8 up to 60° and R/W from 0.1 to 0.5 are presented in Appendix 
A. To duplicate this information for the other three 
meoobems, all that is needed is three additional computer 
Burns. The sample input deck of figure 13 of the text may 
be duplicated and used with a single change for each new 
problem. The entry on card 4 in columns 11-15 will be 
changed to 2, 3 and ¦ for each successive run. 

Each computer run will require 12.5 hours and will 
provide the raw data (a deck of punched cards, each listing 


the problem parameters and corresponding SCF) paralleling 


> 





fee data plotted in Appendix A. A simple plotting routine 
will be needed to plot the data. 

Πα σου ποπ ο the text, we pointed out that the 
SCF at two points are of interest for problem 4, namely the 
Beint at the bottom of the notch and also the opposite point 
che straight boundary. As presented, this program 
calculates only the first of these for each new problem. 
To make the information more complete, an investigator 
Ma tada the card "Y(15) = SSINT(19,1)" following PP 0930. 
Following SS 00090 add the card "SCF2 = SCF * Y(15)/Y(14)"; 
change present card SS 00130 by adding ",SCF2"; change 


present card SS 00140 to "FORMAT (6D15.8)", 
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APPENDIX E 


PROGRAM USER'S MANUAL 


This user's manual is intended to be used with the program 
Sted in Appendix F. 

A user familiar with neither computer coding or the Finite 
Element Method (FEM) of numerical solution can, with the aid 
er this manual: 


П. execute the program listed to solve the four specific 
problems presented in the body of this thesis. 


a er familiar with both computer coding and FEM, in 
addition to the above,can: 
Eebenge the boundary conditions applied to the structure. 
p change the shape of the notch in a minor way. 
4. change the program dimensions to suit the user's needs. 
pee modify the program to solve a different, parameterized 
problem. 
Ш. TO EXECUTE THE PROGRAM 
The job control language cards needed are listed in 
Figure El. No explanation of their make-up will be attempted 
Lis manual other than to comment that some of the cards 
listed beginning with //FTXXF001 provide the direct access 
storage needed by the program. Others provide linkages 
i ed to print, plot, and punch cards. 
A sample input deck is included in Figure El. The entries 
L he cards are listed in Table El in detail. Figure Ee 
(identical with Figure 9 in the text) is a sketch of the 


structure considered, 
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Fig. E2. Structure Faces 


Figure E3 (identical with Figure 8 in the text) shows the 
four problems solved. The user must make some selections 
and can change other. numbers; however some entries on these 


cards should not be changed. They are indicated by: 


User must choose - -( )- 
User may choose ( ) 
User may not choose «ү )* 


Some additional explanatory material may be found in the 


thesis by Adamek [1]. 


moves to Accompany Figure El and Table El. 
Card 4: if NPLOT= 1 or 2, only a single problem is solved. 


ΓΙ. ἱΠΠ το avola unnecessarily drawing the same mesh 
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NTN ONLY MIGOT moczsrleatıon 225 times. This also 
provides the capability to stack as many problems as 
desired in one data run. If NPLOT Z 0, at the end of 
πη’ Τοῦ елдо problem, control returns to the 
reading section of MAIN to begin a new problen. 

Program PLIMEG has an ordering of SEL's with which the 
present use conforms. This means that here we use only 
SEL's numbered 1, 2, b, and 6. Cards 5, 6, 7, and 8 
ТОКОЛ О SEL CS numbered 1, 2, 5, and 6, in that order. 
Cards 5-8: Element size increases with coordinate 
direction for positive EXPX and EXPY. For the smallest 
elemento lle ar the base of the notch EXPX should be 
positive and EXPY should be negative. If the values 
shown ere not used the stresses calculated may be farther 
from the converged result than 1% shown in the body of 
5 thesis. 

Cards 5-8: ROWS must be the same for SEL's 1 and 5 and 
ου ELIS 2 and 6. Likewise for EXPY. 

Cards 5-8. CLMNS must be the same for SEL's 1 and 2 and 
pees and 6. Likewise for EXPX. 

Cards 5-8: Dimension statements in the program limit 
Τα ΠΡ α5 follows; ROWS: SEL1 + SEL2 < 9; 
eu: obi + SEL5 < 15. 

Pd card cOlUuMmis not listed are not used and should 
per left blank. 

ΙΙ; For a discussion of NGP see reference [16], 


Chapter 8. For this formulation NGP = 2 is adequate. 
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NEE SET 


Πο η πος ENFUT DECK 


















Card Variable 
Number Name Explanation 

problem parameters 

(RW) radius/gross width ratio 

(DW) net/gross width ratio 

EIA) Notch semi-angle (degrees) 
identification 

-(TITLE)- | as desired by user 
additional identification 

-(TITLE)- | as desired by user 


mesh parameters 






*(NSEL)* number of SEL's 






* (NPT) * number nodes per element 





(ШЕШ ОН У bir qicator for problem 1 









2m Кл шот tor problem e 


3 indicator for problem 3 









4 indicator for problem 4 





=(NELOLT)= 1) O indicator for no plots 









Па сао or printer plot 












μαι ο τος CALCOM PLOT 





super element deck 


1-5 *CSELNR)* | SEL identification number 

6-10 (ROWS ) number of columns (divisions 
in y direction) in this SEL 

11-15 (CLMNS) number of columns (divisions 


ου τας ο 
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бага cara 


Number 





16-20 





21-25 


26-30 





51255 
36-40 





41-45 





46-50 





51-60 





61-70 





11-15 


26-30 


ον ο 






Columns 











TABLE El 


Variable 
Name 


* (Node 
* (Node 


* (Node 





* (Node 


*(TYPE)* 






*( CON) * 


(EXPX) 


09429 


* (NMAT ) * 


(NSTRES ) 







ДО 





B)* 





m 


Ву 













* (МРТВ) * 












Continued) 


Explanation 
upper left corner node number 
lower left corner node number 
lower right corner node number 
upper right corner node number 
Ша е аЛ делсе ication number 


number of boundary nodes in this 
SEL 


code for’ interconnection between 
SEL — for use see thesis by 
Adamek [H] 


exponential coefficient 
xX direction 


exponential coefficient 
τι ο. ος 96η 


FEM parameters 


number of materials in the 
stmueture 


2 dimension assumption made: 
0 for plane stress 
T for plane strain 


number of Gauss points used in 
the integration 


TOS 





Computer time required may be estimated by knowing that 
the compile and link steps require approximately two minutes 
total and each problem executes in three minutes. Core 
requirement with this configuration is 690K. The core 
requirement may be reduced to 615K if cards with serial 
numbers GG 00530, GG 00740 and GG 00750 are removed from the 
deck. This removes the plotting capabilities from the 
peogram. If NPLOT= lor 2 is used, the program termimates 
mee end of file on unit 5" (the card reader is empty) after 


the last problem is solved. 


2. TO CHANGE BOUNDARY CONDITIONS 

Computations for the load boundary conditicns and codes 
indicating the displacement specifications are entirely 
memuained in subroutine COORD. The serial numbers on ail 
FORTRAN statements used to calculate boundary conditions 
begin with the three digits EEX or EEY. These instructions 
assume that the user can calculate externally applied forces 
B and Ey to describe the loading desired. 

Nodal displacement specifications may be imposed in any 
of several ways. Each specification is indicated by a code 
assigned to one specific node. The array NBC (454,2) is 
Bed to transmit the coded information. NBC (1,1) is the 


u dio placement specification 


node number at which the I 
lies anda NBC (1,2) is the code for that specification 


Faken from the following list. 
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1: A roller on an x axis (i.e.,v 0) 


КОШ roller On a y axis (i.e., u 0) 
3: <A fixed node (i.e., u = v = O) 

A cumulative tally of the number of nodes where a 
displacement specification has been imposed IO aimed 
in the variable NPBC and the tally appears in the output. 

At least one nodal displacement must be specified in 
n Of Che bwo coordimate «directions іп order to prevent 
rigid body motion which might result in meaningless outputs. 


The cards which must be modified are serialized EEY. 


The actual changes to be made are left to the user. 


fee LO CHANGE THE SHAPE OF THE NOTCH 

Only a minor change can be accomplished easily. The 
EUH shape is modified but no modification of the SEL 
arrangement is contemplated. That is, we only wish to 
relocate SEL nodes 1 through 7 (Fig. E4, identical with 
Fig. 12 in the text) and subsequently to generate an 
ШО priate FEM mesh. It is assumed that the same problem 
parameters (R/W, D/W, 8) are used. 

I uubproubiHe PROBL, tche x and y coordinates of nodes 
l through 7 must be re-coded. The FORTRAN statements that 
must be changed can easily be identified by the user. Some 
computational short cuts may be useful under certain 
Ertlons. Before introducing the shortcuts, however, the 
Meer must understand the data string output by this 


Bubroutine. 


1.05 
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We refer to each data string written on direct access 
device 10 as a node card which has the format (315, 2D25.16). 
The first three numbers (integers) are the number of the SEL 
εδρα) Me node number (variable K) and the short- 
Entting code (Variable I or L). The other two numbers 
(Real*8) are the x and y coordinates of node K (variables 
X(K) and Y(K) respectively). The node cards are grouped by 
SEL number and the groups arranged in increasing numerical 
order. Within each SEL group the cards are arranged 
ОООО Ор vo the SEL connectivity. That is, the node in 
the upper right hand corner first, followed by one card for 
each node used to describe the SEL, ordered in counter- 

E ciwise fashion. In this program SEL's 1, 5, and 6 have 
Meenodes each; SEL 2 has 4 nodes. 

ШОШО О ШЕЛ О ОШ seme familiarity with PLIMEG and PLISOP 
is assumed. The coded numbers are described and/or 
merculated as shown in Table De. 

The shortcutting code referred to above may be used if 
the SEL boundary between any two corner nodes is a straight 
line. L=5 for 12 node elements (or L=4 for 8 node elements; 
not used in this program) causes PLIMEG (the mesh generator 
itself) to calculate the coordinates of these intermediate 
modes. Therefore the coordinates need not be calculated in 
Subroutine PROBLM. Card serial BB 00790 is an example of the 
use of the shortcutting code. Note that where the short- 
cine code L 1S not used, I= 0 must be specified as has 


been done in card serial BB 00830. 
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The structure boundary nodes along the circular arc 
(points 1 through 4) are corrected (see page 38 of text) 
in ο... uae COORD. A user modifying Subroutine PROBLM 
must also recode the corrections made by cards serialized 
BEZ. The variables CORD(N, 1) and CORD(N, 2) are the x and 


y coordinates of node N respectively. 


4. TO CHANGE DIMENSIONING 

ο ο. r oli oC Oo am Table δὲ include all of the 
controllable dimension statements in the program. Each is 
exact reproduction of its parent card іп the program 
except that the dimension numbers, where appropriate, have 
Meen replaced by a code: NR2 through NR15. By calculating 
the numbers represented by these codes, preparaing a new 
dimension deck and replacing the parent cards in the program 
with the new cards (on the basis of serial number) a user 
may tailor the program size to meet his needs. 

The designations included in the name column are variable 
names used in PLIMEG and PLISOP. If several different 
Persions of the dimension deck might be required, the computer 
terminal can Be advantageously used. This can be done by 
mains the coded deck into the system, then using up to 
fourteen change commands to change the codes to the previously 
calculated numbers, and then ordering a new punched deck. 

If several problems requiring different dimensions are 
to be solved, they can all be solved using a program with 


mee largest of the dimensions. From the standpoint of 
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efficient computer usage, however, it is better to dimension 
Che program to fit each job. 
A rough estimate of core storage required (in K BYTES) 


may be arrived at by making the following computation: 


2 * NOH * NO10 


12 + 180 


CORE = 


5. TO MAKE MAJOR MODIFICATION 

A major modification to the program involves, primarily, 
a complete rewrite of subroutine PROBLM. The guidance of 
Section 3, together with a good working knowledge of PLIMEG, 
сап be extended to accomplish this. The user may be guided 


by carefully studying the preceding portions of this appendix. 
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MAXMAT 


MAXNJT 


NGCOL 


NGROW 


MAXNEL 


NJT 


NGP 


NBAND 


TABLE El (Continued) 
EXPLANATIONS OF 
DIMENSIONING CODES 


The maximum number of materials to be 
considered 


The maximum number of nodes to be considered 


The maximum number of columns in the basis 
Ες eid 


The maximum number of rows in the basic 
Er 


The maximum number of elements to be 
considered 


The maximum number of nodes per element 
π᾿... “πι; 


The maximum number of Gauss points to be used 
in the integrations (2 through 5) 


NGROW * NGCOL * NO 11 
(NGCOL + 1) * NGROW 
2 * MAXNJT 


The half bandwidth of the stiffness matrix 
calculated from: 


IF NO11= 4, NBAND= 2 * (ROWS + 3) 
IF NO11= 8, NBAND= 2 * (3 * ROWS + 5) 
IF NO11= 12, NBAND= 2 * (5 * ROWS + 7) 


Note? ROWS 15 the maximum number of rows of 


elements in the completed mesh 
MAXNJT + 1 
2 * NO11 


С Oa 
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COMPUTER PROGRAM LISTING 


THIS LINE 


RECORD OF SUBSEQUENT MODIFICATIONS SECULCD FOLLY 


END CE RECORD OF MODIFICATIONS. 
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